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Abstract 

The self-consistency equations for the independent order parame- 
ters as weh as the free energy expression for the mean-field RVB model 
of the spin- 1/2 Heisenberg Hamiltonian on the anisotropic triangu- 
lar lattice is considered in the quasi-one-dimensional approximation. 
The solutions of the self-consistency equations in the zero-temperature 
limit are in fair agreement with the previous numerical analysis of the 
same model by other authors. In particular, the transition from the 
ungapped ID-RVB state to the gapped 2D-RVB state occurs at an 
arbitrarily weak transversal exchange (J2 — )• 0) although the amount 

of the gap is exponentially small: exp 1 where Ji is the 

longitudinal exchange parameter. The structural consequences of the 
formation of the 2D-RVB state are formulated by extending the fa- 
mous bond order vs. bond length relation known for polyenes (one- 
dimensional Hubbard chains). Analytical estimates of this effect are 
given. 
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1 Introduction 

The RVB state originally proposed by Pauling [Ij for describing the struc- 
ture of the benzene molecule is being sought in many materials after An- 
derson's p] conjecture that it represents the ground state of cuprate-based 
high-temperature superconductors. The recently obtained [3| CuNCN phase 
whose structure is represented in Fig. [T] had been proposed as a candidate 
for an RVB ground state spin liquid [4J due to frustration of the effective 
exchange in the ab-plane where the 1/2 Cu^"*" local spins form an anisotropic 
triangular lattice. The material had been subject of a series of measure- 
ments of its magnetic susceptibility, electric resistivity, heat capacity, also 
elastic neutron scattering as well as of ESR, NMR relaxation, and muon spin 
resonance (all vs. T) [3] [S] (6]. Although the neutron scattering (complete 
absence of the magnetic signal) as well as the susceptibility measurements 
(temperature independent paramagnetism above 80 K changed to approxi- 
mately activation decay of the susceptibility below this temperature) strongly 
indicate the transition between the ID-RVB regime at higher temperatures 
to the gapped 2D-RVB regime below 80 K the issue remains controversial [6] 
since the NMR and the yuSR data so far better fit into a model of inhomoge- 
neous spin-glass-like ground state. In order to solve this dilemma we try to 
find out in the present paper what might be the structural concequences of 
the formation of two spin-liquid phases which is going to help in obtaining 
direct experimental evidence of this picture. 

2 RVB mean-field analysis of an anisotropic tri- 
angular lattice system 

2.1 Hamiltonian 

A close inspection of the structure (Fig. [1]) reveals that each Cu^"*" ion can 
be effectively antiferromagnetically coupled to two of its neighbors forming 
a chain while somewhat weaker antiferromagnetic coupling with four more 
neighbors from two adjacent parallel chains results in a Heisenberg model on 
an anisotropic triangular lattice with the Hamiltonian: 
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Figure 1: A look into the CuNCN crystal structure. A stronger Ji extends 
in the a direction; a somewhat weaker J2 extends along the 6 ± a directions. 
The weakest J3 extends in the c direction and is not considered in the present 
paper. 
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where the couphng vectors r take three values Tf, i = 1 3; Ti = (1, 0); T2 = 
(|,^); T3 = (|,— with the interaction of the strength Ji along the lattice 
vector Ti (two neighbors) and with a somewhat smaller strength J2 along the 
lattice vectors T2 and T3 (two neighbors along each). This is precisely the 
setting for which Yunoki and Sorella ^ proposed that two (different) spin- 
singlet RVB (s-RVB) states are formed at different temperatures depending 
on the amount of anisotropy ^ on the basis of their VMC calculations. Later 
Hayashi and Ogata [8] reproduced this result within a mean-field treatment 
which we basically follow here. 



2.2 Equations of motion and self consistency equations 

Hayashi and Ogata [8J base their analysis of the Hamiltonian eq. ([T]) on 
returning to the electron representation from the spin representation by the 
standard formulae: ^ 

Si = ^cl^(^al3Cip, (2) 

where c^(Qo-) are the electron creation (annihilation) operators subject to the 
Fermi anticommutation relations; cr^^ are the elements of the Pauli matrices 
and the summation over repeating indices is assumed. For the latter one can 
derive equations of motion based on the Heisenberg representation in which 
each operator obeys the following equation of motion: 

ihA = [A, H] (3) 

where [, ] stands for the commutator of the operators and the dot-on-top 
symbol for the time derivative. Applying this to the creation and annihilation 
operators c^^{cra) and performing commutation, mean-field decoupling and 
Fourier transformation as done in j3l E] results in mean-field equations of 
motion for these operators : 

ihci,^ = -- ^ Jt^t cos(kr)cka -| Y.r J-r^T cos(kr)c+,^_^ 

T 

^^c^. = ^ 5Z ^°^(^^)^k/3 +|E.^tA; cos(kr)c_k-<x (4) 

r 

These reduce to the set of 2 x 2 eigenvalue problems for each wave vector k: 



Ak -a / V / V 
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with 

Ck = -3 ^T^T cos(kr) 

Ak = 3 ^ cos(kr) (6) 

T 

(summation over r extends to ±rj; i = 1^3) which resuhs in the eigenvalues 
(quasiparticle spectrum) of the form: 

= (7) 

whose eigenvectors are combinations of the destruction and creation oper- 
ators with the above Bogohubov transformation coefficients u^, v^. These 
equations resuh in the self-consistency equations of the form: 




for six order parameters , A,- defined as 

At = (CraCr+r/?) • (9) 

It is remarkable and important for the subsequent treatment that the order 
parameters C,t are in fact bond orders for the corresponding pairs of atoms. 

2.3 Free energy 

Following Ref. \9\ one can write immediately the free energy in terms of the 
above order parameters: 

m - ~ E (2 cosh ) + I E Ml 4^Jr l^rl' (10) 

where 6 = ksT and summation over r extends to ±rj;i = 1^3. Minima 
of this expression with respect to C,t and A^. correspond to various possible 
states of the system. 
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3 Simplified RVB on the anisotropic triangular 
lattice 



The numerical analysis [8] shows that in agreement with general theorems 
[lOj the order parameters satisfy additional phase relations 

arg - arg A^.^ = ±| (11) 

which allows to reduce the number of order parameters to only two [TT| H]: 

V2^ = ^^,;V27]=\A^,\ = \A^,\ (12) 

the first responsible for establishing the gapless ID-RVB state within the 
chains and the second one for opening a gap and establishing the 2D- RVB in 
the transversal direction and other three order parameters set to be zero. It 
can be shown that this particular choice of the phases leads to the following 
form of the fermion quasiparticle spectrum: 

El = 18 (^Ji^^2 cos^iK) + JW + JW cos(k^) cos(^y3)) . (13) 

This form of the quasiparticle spectrum allows for a simple analysis. Obvi- 
ously in the ID-RVB state {rj = 0) the spectrum has a vanishing dispersion 
in the y-direction: 



=3^2 |J^cos(k^).| (14) 

It is gapless along the lines k^, = ±| and contains ridges at k^. = 0,±7r 
extended in the y-direction. The corresponding density of the quasiparticle 
states (qDOS) is depicted in Fig. 2 (red). The divergence on the qDOS: 



on the upper border of the spectrum is obviously due to the abovemen- 
tioned ridges in the spectrum of the quasiparticles. On the other hand the 
low-energy behavior of the system in the ID-RVB state is controlled by a 
constant qDOS at the zero energy, which is perfectly reflected in the temper- 
ature independent paramagnetic susceptibility in the respective temperature 
region. 
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qDoS 




Figure 2: The schematic densities of the fermion quasiparticle states in the 
ID-RVB (red) and 2D-RVB states (blue) of a system of 1/2 Heisenberg spins 
on the anisotropic triangular lattice. The 2D-RVB qDOS is built for the 
effective anisotropy parameter a = | (see the text for details). In addition 
to the gap opening one may observe a slight shift of the energy at which the 
upper singularity occurs. 
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Table 1: Critical points of the quasiparticle spectrum in the gapped spin 
liquid 2D-RVB state of a system of 1/2 Heisenberg spins on the anisotropic 
triangular lattice. (ka;,kj,) stands for the coordinates of the critical point in 
the Brillouin zone, n is the degeneracy - total number of points of the given 
type; other entries are self-explanatory. 



n 



point type 



-2,0,2 



18 {jfe + 2J2V) 
18 (ite + 2j-v^) 



maximum 



maximum 



2t2 



saddle 



(±7r,0) 



2 £-2 



18J/^ 



saddle 



(±E ±^I—) 



saddle 



2^2 



saddle 



f± arccos 



■2J252 



±^1 



2^2 



1 I 3 
^ + life 



saddle 



f± arccos 



2W 



±— 1 



2^2 



4 J?F 



minimum 



f± arccos 



.0) 



1 72„2 

/ 4 



minimum 



The transition to the 2D-RVB state is accompanied by a significant re- 
organization of the quasiparticle spectrum. Both the ridges at k^. = 0, ±7r 
and the degeneration lines at k^; = ±| disappear and are replaced by the 
critical points whose characteristics are given in Table [T] One can realize 
that in the 2D case the qDOS is mainly contributed by the logarithmic van 
Hove singularities due to the saddle points of the spectrum on top of the con- 
stant contributions coming from the minima and maxima of the quasiparticle 
spectrum. Physically it must be important that the values of the energy cor- 
responding to the minima and maxima differ from the saddle point energies 
only in the second order in a small parameter of efficient anisotropy: a = ^ 
so that one ultimately cannot expect anything, but some widening of the 
logarithmic peak of the qDOS on the upper bound of the spectrum. The 
most important changes occur at the lower bound of the spectrum, where a 
gap opens. As one can see from Table [T] the singular spectral weight must 
concentrate at the energy of the lower saddle points: 8-^/2^/2^, although the 
lower boundary of the spectrum is smaller than this value in the fourth order 
with respect to efficient anisotropy. 
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On the other hand we notice that the wave vectors of the critical points 
of the quasiparticle spectrum in the vicinity of its lower boundary either 
have the k^^ component equal to ±^ or differ from these values in a higher 
order of efficient anisotropy. Following [TT| we notice that in this range of 
wave vectors the dispersion of quasiparticles in the y-direction is negligeably 
small. This brings us to the idea that one can hope that neglecting the y- 
dispersion when calculating the integral characteristics of the system does 
not affect the precision catastrophically. This may be considered as a quasi- 
one-dimensional approximation for the spectrum, which then takes the form: 



= ^18(J2ecos2(k,) + J2V). (15) 

Inserting this in the standard definition of the density of states we obtain as 
expected: 

9{^) = (16) 

which corresponds to the quasiparticle band ranging on the energy scale from 



3-\/2^2^ to 3-\/2a/ Jli"^ + J^V^ with its lower boundary being as explained 
above somewhat higher than the lower boundary of the exact spectrum and 
with the upper boundary being located between the upper boundary and 
logarithmic peak of the exact spectrum. 

With use of the qDOS eq. ( fT6l) one can easily write the explicit expression 
for the free energy as relying on the general expression eq. ffTOj) . It reads as 
follows: 

F = 6 Jie' + 12 Jar/' -26 j g{s) In (2 cosh (^^) ) de. (17) 

We do not expect that our results obtained in |1] with use of the high- 
temperature expansion: 

E]iW , . ^ f iZ-Ek^"^ 



In 2 cosh { ^ ] ] ^\n2+ , ^ , , ^ 

V \29 J J 2\2d J I2\2d 

which allowed to perform integration over the entire BZ without neglecting 
the y-dispersion are going to change. By contrast, in the low temperature 
regime we first rewrite 
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^ In ( 2 cosh ( ^ 



= £ + + exp(-i)) 

and immediately obtain the ground state energy for the 2D-RVB state: 



F{9 = 0) = 6Ji^^ + 127277^ - / g{e)ede 



The integral is done analytically |12] : 

6\/2 



TT 



72^2 + J2VE(fc) 



where E(/c) is the complete elliptic integral of the second kind of the modulus 
k given by: 

72f2 

This result is not unexpected since it has a form characteristic for one- 
dimensional systems [i3\. Taking derivatives with respect to the order pa- 
rameters and setting them equal to zero results in self-consistency conditions: 

1 = -— ^ K (fc) 

which are remarkably similar to the self-consistency conditions in the 
one-dimensional Hubbard problem, the first being one for the bond order 
and the second being analogous to that for the gap or magnetization with 
the J2 parameter taking part of the interaction parameter U of the Hubbard 
model and 3a/2Ji^ being the effective one-dimensional bandwidth. In the 
ID-RVB state the first of the conditions eq. (ITSl) yields the amplitude of the 
order parameter ^ reached at the zero temperature: 

which is in perfect agreement with the numerical result of [8]. Inserting this 
in the second of the two conditions eq. ( fTSj) . neglecting the terms with con- 
taining T] as compared to those with ^ in the sums, and using the logarithmic 
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asymptotic form of the complete elliptic integral of the first kind K we arrive 
to the self-consistent field-like estimate of the for t] and for the gap in the 
2D-RVB state: 



2v^Ji / 2Ji 



3V2J2V 



J exp 

12Ji / 2Ji 
exp — 

TT V J2 



at zero temperatue. The latter result is in a fair agreement with the numerical 
study pn] where the pre-exponential factor in the gap was estimated to be 
3.50 as compared to 12 /n ^ 3.82. Although the factor in the exponent was 
estimated to be 1.61 in [11] against our estimate of two, the general form 
of the dependence of the characteristics of the model on its parameters is 
reproduced. These results in agreement with the numerical results [TT] and 
general behavior of one-dimensional models with interaction and show that 
at the zero temperature some nonvanishing value of the 1] order parameter 
and the energy gap appear at arbitrarily weak interaction J2 so that no 
critical point with respect to the anisotropy J2/J1 should be expected at 
zero temperature. 

A further move consists in inserting the above expression for 7] in the 
logarithm and retaining the terms proportional to Afc'^; A = In p; fc'^ = 1 — /c^ 
in the equation for ^ so that K{k) — D (k) ^ 1 — Afc'^/2. By doing so and 
retaining the terms up to second order in 7] we obtain: 

which represents the estimate of the bond order variation in the 2D-RVB 
state as compared to the ID-RVB state. 

This finding is in a fair and remarkable agreement with the numerical 
result of [11] where it was shown that in the region where the 2D-RVB state 
develops {1] 7^ 0) the C, parameter manifests a very weak depletion as com- 
pared to its ID-RVB [rj = 0) value. Despite the fact that it must be not 
particularly strong this depletion can manifest itself in a geometry change, 
which in principle could be observed. The subsequent reasoning following 
the lines proposed in [15] shows how it can tentatively look like. 

Assume that the lattice contribution to the total energy per copper site 
can be harmonically approximated as a function of the separation p between 
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the Cu atoms in the a-direction (the shortest such separation in the struc- 
ture): K{p — Pq)^/2, where K is an effective elastic constant and pg is an 
equihbrium separation to be observed if the spin contribution eqs. f fTOj) . f fT7]) 
to the energy are turned off. The spin contribution to the lowest order comes 
from the "kinetic energy" term and is given by — (6/2)/7tJi^. Assuming the 
dependence of the effective exchange integral on the interatomic separation 
in the formJi = Jiq + J((p — p^) we easily arrive to the analog of the famous 
bond order vs. bond length relation [16j for the RVB states: 

p-Po = ^{J[0/K, (21) 

where p is now the equilibrium interatomic separation in the presence of the 
spin contribution to the energy. Obviously, the separation is going to change 
according to the sign of the derivative of the effective exchange integral with 
respect to the interatomic separation increase {J[ < 0). The latter condi- 
tion is, however, natural, in the frame of the standard conception of the 
sources of the antiferromagnetic exchange. Indeed, these appear as a result 
of perturbative treatment of the one-electron hopping in a strongly interact- 
ing regime. Then one has: Jj = {Atf)/U where is the intersite one-electron 
hopping parameter along the respective hopping vector Tj and U is the on- 
site electron-electron repulsion parameter. Assuming a linear dependence of 
the hopping parameter on the interatomic separation: t = + t' {p — p^) 
we arrive to an estimate J[ = {8tot')/U < since one can easily see that 
the multipliers to and t' must have opposite signs. Combining eqs. f l20|) and 
(12T]) we arrive to the estimate for the variation of the equilibrium interatomic 
separation in the 2D-RVB state: 

8p = -6J2J[vVJiK 

which immediately shows that the lattice parameter a in CuNCN must man- 
ifest the same trend as the 2D- RVB gap {J[ < 0), although with somewhat 
damped amplitude due to square in a small quantity i]. 

4 Conclusion 

In the present paper we succeeded in obtaining analytical estimates for the 
parameters (order parameters, energy gap) of the 2D-RVB state and possible 
structural concequences of the variation of these parameters under a ID to 
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2D-RVB transition conjectured recently [B] to be responsible for the observed 
temperature behavior of the magnetic susceptibility in CuNCN. 
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